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Abstract 



We use Clebsch potentials and an action principle to derive a closed sys- 
tem of gauge invariant equations for sound superposed on a general back- 
ground flow. Our system reduces to the Unruh (1981) and Pierce (1990) 
wave equations when the flow is irrotational, or slowly varying. We illustrate 
our formalism by applying it to waves propagating in a uniformly rotating 

fluid where the sound modes hybridize with inertial waves. 
PACS numbers: 43.20B, 43.28Py 
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Typeset using REVT^ 



I. INTRODUCTION 



Studies of sound in a stationary fluid usually start from from the wave equation whose 
derivation appears in all elementary texts j^. When the fluid is moving, however, finding 
an appropriate generalization of the wave equation is not always possible. An exception is 
the special circumstance that the equation of state is barotropic and the background flow 
irrotational, but not necessarily steady. In this case a particularly attractive equation was 
derived by Unruh . Unruh's equation was later rediscovered and further popularized by 
Visser It coincides with the equation obeyed by a relativistic scalar field propagating 

in curved space-time. The space-time geometry is governed by the acoustic metric which 
depends on the background flow velocity and on the local fluid density and speed of sound. 

The curved space-time interpretation of the wave equation is rather more than a mathe- 
matical curiosity. As well providing an attractive analogy with some aspects of general rela- 
tivity 0, one can use the geometric formalism for ray tracing, and to produce a straightfor- 
ward and systematic derivation of various conservation laws associated with acoustic energy 
and momentum 

Unfortunately most fluid motions occurring in nature are not irrotational. It is therefore 
desirable to explore the possibility of extending the acoustic metric equation to a wider class 
of flows. A valuable step in this direction was taken by Pierce |^ who, without assuming that 
the background flow was irrotational, derived an equation which is in appearance equivalent 
to the acoustic metric equation. He did, however, assume that the background flow varied 
slowly over the length and time scale of the sound wave. As a result of his approximations, 
his dependent variable is not quite the velocity potential appearing in the acoustic metric 
equation. It gives the fluctuating velocity only up to a correction whose magnitude depends 
on the space and time inhomogeneities of the background flow. Pierce did not attempt 
to characterize the correction beyond estimating its size and showing that it was small 
in most regimes of interest in acoustics. He did show, however, that in special cases his 
equation reduces to known exact wave equations. In particular, for steady irrotational flow 
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it reduces to Blokhintsev's equation |TyjTT|, which is a special case of Unruh's. This and 



other features led Pierce to conjecture that his equation is of wider applicability than his 
derivation suggests. The purpose of the present paper is to show that this conjecture is 
correct. 

We use Clebsch potentials and an action principle to derive the equations of motion for 
small perturbations about a general barotropic flow. Our principal result is an exact and 
concise expression for the small correction to potential flow, and a simple equation of motion 
obeyed by it. This, coupled with the Pierce equation, provides a closed system for wave 
propagation in a general inhomogeneous and unsteady background flow. The condition that 
the correction to potential flow be ignorable, and thus the acoustic metric equation accurate, 
is that the frequency of the sound be appreciably higher than the local vorticity. There is 
no requirement that the spatial inhomogeneity be small. Even if the frequency condition is 
violated, we can still study the wave motion, but with a more complicated system of partial 
differential equations. 

The paper is organized as follows: In section |T| we review the action principle for the 
Clebsch formulation of barotropic fluid mechanics. In section ^ we discuss the infinite 
family of conserved quantities that generate global gauge transformations on the Clebsch 



potentials. These are needed in section |I^ where we consider first order perturbations to 
a background flow, and identify two gauge invariant combinations of the potentials which 
have physical significance. Sections and ^ derive a closed system of equations for these 



combinations. Finally section |VI1| illustrates our formalism by applying it to waves in a 
uniformly rotating fluid. 

II. EULER EQUATION 

In this section we will review the Clebsch potential approach to fluid dynamics. The 
Clebsch formalism has the advantage that the equations of motion may be derived from an 
action principle |T2[ , and with an action principle conservation laws are related to symmetries 
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by Noether's theorem. 
We begin with 

S = Jd'x j-lpv^ - 0(p + V • (pv)) + p/?(7 + (v ■ V)7) + «(p)} . (2.1) 

Here p is the fluid mass-density, v the velocity, and u{p) the internal energy density. This 
is the customary expression giving rise to irrotational fluid dynamics (see, for example [0) 
— but with an additional term containing new fields: f3 and 7. The variable f3 may be 
thought of as a Lagrange multiplier enforcing the Lin constraint |13| that there be a label 
(7) painted on the particles permitting us to distinguish one from another^. 
Requiring that S be stationary when we vary v gives 

- pv + pV0 + P/3V7 = 0, (2.2) 



or 



v = V(/> + /?V7. (2.3) 



This is the Clebsch representation |l5| , p^ of the velocity field. It allows for flows with 
non-zero vorticity 

= V X V = V/9 X V7. (2.4) 

We use ( p.3|) to algebraically eliminate the v in in favour of the Clebsch potentials 0, 
f3, 7. This leads to a new action |jl2 



5new = / rfV |^p(V0 + /5V7)' + P(0 + Pi) + m(p)| . (2.5) 

Varying the remaining variables in ( p^.5| ) gives the equations of motion 



^Lin originally employed three Lagrange multipliers /3i,2,3 leading to the conservation of three 
Lagrange co-ordinates, 71231 which served to label the material particles uniquely. As shown by 



Seligar and Whitham [14|, only one is really necessary. 



(50 : p + V ■ (pv) = 0, 

6(3 : p(7 + (v • V)7) = ^ 7 + (v ■ V)7 = 0, 

^7: 9t(p/?) + V(vp/?) = ^ /5 + (v ■ V)/? = 0, 

5p: it;2 + + /?7 + /i = O, (2.6) 

where, in the last hne, p = du/dp is the specific enthalpy. We see that the values of both (3 
and 7 are advected with the motion. 

We now show that these equations imply the Euler equation for the fluid. Apply Vj 
to the last line in ( p.6|) and add and subtract /3Vi7, so that the second term is the time 
derivative of the velocity 

VkViVk + ^t{V^(t) + /3V,7) - /^Vi7 + 7V*/? = -V^p. (2.7) 

In other words 

dtVi + WfcVit;fe - /5V,7 + 7Vi/5 = -Vip. (2.8) 

The second, third, and fourth terms on the left hand side now need to be taken care of. 
Write 

Vk^iVk = VkVkVi + Vk{ViVk - VkVi), 

= {v-V)vi + Vk{Vi(3Vkl -Vk(3Va), 

= (vV)vi-^ViP + pVa, (2.9) 

where, in the last line, we have used the convective constancy of (3, 7. 
Inserting (|2.9| ) into (p.8|) we find 

v + (v- V)v = -Vp, (2.10) 

which is Euler's equation. 
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III. GAUGE TRANSFORMATIONS 



In three dimensions any vector field may be locally represented in the form [jT5|,|T6 



v = V0 + /?V7. (3.1) 

Given a velocity field v, however, the potentials 0, (3 and 7, are not uniquely determined. 
This indeterminacy is usually described as being due to a gauge invariance, but it is more 
analogous to the residual freedom to make time independent gauge transformations that 
survives after we fix the Aq = gauge in Maxwell electrodynamics. As in that example, 
once we have a made a choice of the potentials, 0, (3, and 7, at any particular time, their 
subsequent evolution is uniquely determined by the equations of motion ( p. 61) . 
We can relate the gauge invariance to conservation laws. From 

p + V ■ pv = 0, 
/3+(v-V)/3 = 0, 

7 + (v • V)7 = 0, (3.2) 

we deduce that 

' pr{f3,-i)d^x (3.3) 



is independent of time. Here is an arbitrary function of the variables j3 and 7 with position 
independent coefficients. 

Now any action that contains only first order time derivatives defines a Poisson bracket 
and canonical structure. For two functionals Fi 2 of the fields p, 0, (3, 7 at time t we define 
the Poisson bracket {^1,^2} as 

dF2 



dt 



(3.4) 



where the subscript, Fi, on the derivative indicates that time evolution of the variables p, 
0, /3, 7 is derived by varying the action 



S[F^] ^lp{<j> + /37) d^'x - I F^{p, 0, /3, ^)dt. (3.5) 

Such a Poisson bracket automatically satisfies all the usual properties, including skew sym- 
metry and the Jacobi identity 

{Fi, {F2, F3}} + {F2, {F3, Fi}} + {F3, {Fi, F2}} = 0. (3.6) 

In the present case the bracket becomes 

r ,3 /I 5Fi 5F2 SF, SF^ (5 5F, SF^ _ ^ \ 

{F„ F,}-jdx i^-^^ - - p w " ^ " r ^ ^ ^ 

and (p, </)), and (p/9, 7) constitute two canonically conjugate pairs, i.e. 

{p(x),0(x')} = <5'^(x-xO, 

{p/3(x),7(x')} = <5^(x-x'). (3.8) 

We now consider the conserved charge F as the generator of an infinitesimal symmetry 
by setting 

5cl>^{F,cl>}^T-P—. (3.9) 

Similarly 

07 

(3.10) 

The field p is unaltered. This is because F does not contain (f). 

These variations generate an infinite dimensional global (rigid) symmetry group. It is 
a global symmetry because the parameters in JF are required to be independent of x and 
t. The transformations are the extension to Clebsch potentials of the global U{1) phase 
symmetry — > + const, which appears in potential flow, v = V0, where it is generated 
by the conserved charge Q = J p d?x. 

The symmetry transformations leave the Hamiltonian 



H = J {^p(V0 + /?V7)' + d'x (3.11) 

invariant because {F, H} = —{H, F} = ^ = 0. In addition to Poisson-commuting with the 
hamihonian, the conserved charge F generates variations that preserve v itself: 

= 0. (3.12) 
They also preserve the kinetic term: 

which vanishes provided F does not explicitly depend on time. 

It is easy to show that the symmetry group is that of orientation and area preserving 
diffeomorphisms of the 2-plane. It is equivalently the group of non-linear canonical trans- 
formations on a two-dimensional phase space with Darboux co-ordinates (3, 7. Because of 
this we can obtain the finite form of the transformations — as well as confirming that that 
they exhaust all transformations that preserve v — by exploiting the familiar generating 



function methods from classical mechanics [|T^]. 
Suppose that 

d(j) + pdj = d(l) + pdj, (3.14) 

Then 

d{4)-(j)) = pd-f- /3d^, (3.15) 
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and there must exist a VF(7, 7), the generating function, such that 

dW dW 
(t>-(t>^W, ^ = = (3-16) 

Conversely, given a generating function, we can obtain a finite canonical transformation. 
To make contact with the infinitesimal transformations we considered earlier, we let 

^^(3 + $At, 

7 = 7 + 7At, (3.17) 

where ''t" is a notional time parameterizing the change. 
Thus 

dW = (3d-f -{(3 + /3At)d{'y + 7 At), 

^ -At($d-f + (3d^). (3.18) 

Similarly let W = UAt, so that 

dU ^ -$d-f - (3d^, (3.19) 
or, making a Legendre transformation F — U + P'j, 

d{U + (3^)^-$d-f + ^d(3^dF{(3,-f). (3.20) 



In other words 



leading to 



dF dF 



= + 1^ At = 0+ At, 
dF 

P^P-—At, 
dF 

^^-f + —At, (3.22) 



op 



as before. 
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IV. GAUGE-INVARIANT FLUCTUATIONS 



We want to study the evolution of small fluctuations (pi, 0i, 71), superposed on a 
background flow (po, 0o, /^o, 7o)- We will not assume that this background flow is steady, 
only that it satisfies the equations of motion. 

We expand the action ^new out to quadratic order in the fluctuations 

5new = + 5i + ^2 + ■ ■ ■ . (4.1) 

The action 5*1, containing terms linear in the fluctuations, vanishes because of our assumption 
that the zeroth order variables obey the equation of motion. The term quadratic in the 
fluctuations is 

S2 = J d^xi^^povl + pivo ■ vi + pi(0i + /?o7i + A70) + P0A71 + ^"'^ij ' (^-2) 
where vi is shorthand for V0i + PiVjo + PqVji, and 

O d Hi , . 

c' = Po;^, 4.3 

is the square of the local speed of sound. From S2 we can deduce the equations of motion for 
the first-order fluctuating quantities. These equations are not easy to work with, however. 
Because they are advected with the flow, the potentials /3o and 70 which appear as coefficients 
in the equations will generally be time dependent — even if the background flow is steady. 
Furthermore, the overall gauge ambiguity obscures any physical interpretation. It is therefore 
fruitful to seek combinations of the potentials that are gauge invariant and can be expressed 
in terms of physical quantities. For example the ffist order velocity field, 

vi = V01 + /3oV7i + A V70, (4.4) 

is gauge invariant because v is. 
By varying pi in (|4.2| ) we find 



Pi = (01 + vo ■ V01 + /?o(7i + vo • V71) + A (70 + vo ■ V70)) • (4.5) 
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Since 



we can write this as 



where 



/?o + vo ■ V/?o = 0, 

70 + vo ■ V70 = 0, (4.6) 



Po #1 „x 



V'i = 0i+/5o7i, (4.8) 

and 

is the convective derivative. 

The density fluctuation pi, being a physical variable, is gauge invariant. Consequently 
equation ( [4.7| ) suggests that the combination ipi is itself gauge invariant. This is easily 
confirmed: 

^('^i + /?o7i) = + ^71 - - /3o (^^A + ^71 j 

= 0. (4.10) 

We can use ipi to write 

vi=V^i+ei, (4.11) 

where 

ei = AV7o-7iV/?o. (4.12) 

This is a decomposition of the first-order velocity fluctuation into two gauge invariant parts. 
Because sound in a fluid is a scalar excitation, it is natural to identify the scalar field ipi with 
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the acoustic degree of freedom, and ^i, the correction to potential flow induced by angular 
momentum conservation, with a partial hybridization of the sound with other modes. 

Although the vector field .^i has three components, it only represents two degrees of 
freedom. This is because 

^i-Uo = (AV70 - 7iV/3o) • (V/?o X V70) = 0. (4.13) 

Since ^1 is gauge invariant, it should be possible to write it in terms of physical variables. 
In section ^ we will show that it is equal to xi x uq where xi is the particle displacement 
caused by the disturbance. 

V. WAVE EQUATION 

The first-order continuity equation 

+ vo • Vpi + piV • vo + V ■ povi = 0, (5.1) 



dt 

together with the zeroth order continuity equation 



dt 

the equation for pi. 



+ V ■ povo = 0, (5.2) 



Pi = — (5.3) 



c2 dt 

and the decomposition vi = Vipi + ^1, may be combined to give 



If we ignore the .^1, ( |5.4| ) is Pierce's approximate wave equation 



By using (p.2[) again, this can be rewritten as 
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^ + V ■ vo ) ^ ( I: + vo ■ V Wi = V(poV^i), (5.6) 



\dt J c2 \dt 

where each V is acting on everything to its right. 



Ahhough ( |5.5|) may seem more natural, the form ( |5.6| ) has the advantage that it can be 
written a^ 

^ -.d.^g'^-'d^ij^ = 0, (5.7) 



where 

^ Vvo, VqV^ - cn/ 

Equation (|5.7|) has the same form as that of a scalar wave propagating in a gravitational 
field with pseudo-Riemann (Lorentzian) metric g^u. We will refer to g^^ as the acoustic 
metric. The idea of writing the sound wave equation in this way is due to Unruh 

As is customary in general relativity, the symbol g denotes the determinant of the co- 
variant form of the metric, g^^, so detg^'^ = g^^. Taking the determinant of both sides of 
(xH) thus shows that the 4- volume measure y/—g is equal to p^/c. Knowing this, we may 



then invert the matrix g^^^ to find the covariant components of the metric 



" \ vo. 

The associated space-time interval is therefore 

Po 



9,u = ^\ " "1. (5.9) 



ds^ = [c^dt^ - S,j{dx' - vidt){dx^ - vidt)]. (5.10) 

In the geometric acoustics limit, sound propagates along the null geodesies of this metric. 
Metrics of the form ( |5.10|) , although without the overall conformal factor po/c, appear in 



the Arnowitt-Deser-Misner (ADM) formalism of general relativity |T^. There, c and —Vq are 



^We use the convention that greek letters run over four space-time indices 0,1,2,3 with 
while roman indices refer to the three space components. 
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referred to as the lapse function and shift vector respectively. They serve to glue successive 



present case, provided po/c can be regarded as a constant, each 3-space is ordinary fiat 
equipped with the rectangular cartesian metric g\f"''"^'^ = 6ij — but the resultant space-time 
is in general curved, the curvature depending on the degree of inhomogeneity of the mean 

flow Vq. 

This formalism is very pretty, but ( ^.5|) is exact only when the background flow is po- 
tential. Equation (|5.4|) , on the other hand, is valid for a general barotropic flow — but to 
be of use it must be complemented by an equation determining the time evolution of ^i. We 
now derive such an equation. 

We start with the observation that, since P, 7 are convectively conserved, we have 



three-dimensional time slices together to form a four-dimensional space-time [0. In our 



+ (vo ■ V)/3o = 0, 



(5.11) 



at 



and 



d/3: 



+ (vo ■ V)/3i + (vi ■ V)/3o = 0. 



(5.12) 



dt 



Taking the gradient of ( ^.llj ) gives 




(5.13) 



Thus, using the definition ( [4.12|) , 




^ + Vo ■ Vj = -[(vi ■ V)/3o]Vi7o + [(vi ■ V)7o]Vi/5o 

-/3i(Vit;oj)Vj7o + 7i(Vit;oi)Vj/5o 



(-Vj?/'i - ^ij){VjVoi - ViVoj) - (yiVoj)^ij, 



(5.14) 



which can be written as 
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The two equations 



= V^i X a;o - (6 ■ V)vo. (5.15) 



and 



^ = V^i X cuo - (6 ■ V)vo, (5.17) 

form a closed system of equations, containing only gauge-invariant quantities, describing the 
first-order fluctuations about the background mean flow. 



VI. DISPLACEMENT FIELD 

It is not yet clear that, under most circumstances of interest in acoustics, the quantity 
1^1 is a small correction to Vipi. It becomes so, however, once we establish the result 

^1 = xi X uo (6.1) 

where xi is the displacement of a material particle due to the sound wave. Given (|6.1|), we 
see that the order of magnitude of is that of the product of the displacement amplitude 
with the background flow rotation frequency. The fluctuating velocity associated with the 
acoustic field is, on the other hand, of the order of the displacement amplitude times the 
frequency, Q, of the sound wave. ThusQ is smaller than Vipi by a factor of Ic^ol/fi. 

To establish ( |6.1| ) we first define xi(x, t) to be the eulerian displacement field (i.e the 
material point which in the unperturbed reference flow was at time t located at x IS, as a 



Our argument tacitly assumes that xi remains small and oscillating. This is certainly what we 
expect for a sound wave, but, in the absence of viscous damping, many flows with vorticity will 
be unstable to the onset of turbulence, and if the sound triggers such an instability xi will grow 
without bound. 
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result of the perturbation, now to be found at x + xi). We remember that the numerical 
values of the potentials /3, 7, are painted on the material particles, and so move with the 
flow under both time evolution and the creation of an initial perturbation by means of an 
external potential body force. Interpreting this statement mathematically leads to 

Xl ■ V/?o + A = 0, 

Xl ■ V70 + 71 = 0. (6.2) 

From this we may write 

AV70 - 7iV/3o = (xi ■ V7o)V/?o - (xi • V/3o) V70, 
= Xl X (V/3o X V70), 

= Xl X ujq. (6.3) 

We can use (|6.1| ) to rederive the equation of motion for ^1 and so provide a derivation 
of the wave equation that is independent of the use of Clebsch potentials. In their absence, 
though, the origin of the decomposition of the velocity field into the sum of ^1 = xi x 
and the gradient of the pressure potential, V'l, is a trifle obscure. 

To verify that ( |6.1| ) leads to the equation of motion ( |5.14| ) for ^1 we must first establish 



a connection between Vi and the time derivative of Xi. This requires us to describe the 
perturbation with a little more formality. Consider a family v(x, t. A) of adjacent solutions 
of the full equations of motion. The velocity field v(x, t,0) is that of the unperturbed 
reference flow, and increasing values of A correspond to flows evolving from a one-parameter 
family of initial perturbations. By definition the operations of time evolution and variation 
of A commute. 

The position, x(t. A), of a material particle is given by the solution to the differential 
equation 

x(t)=v(x(t,A),t,A) (6.4) 

with suitable initial conditions. Our first order perturbed fields are, in this language, 
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Differentiating ( |6.4| ) with respect to A, and interpreting the time derivative as a convective 
derivative, gives 

vi = ^ + (vo • V)xi - (xi ■ V)vo. (6.6) 
at 

Now, starting from 

^1 = xi X uq (6.7) 

and the convective derivatives 

(ixi 



dt 
dujQ 



vi + (xi ■ V)vo (6.8) 
{uq ■ V)vo - (V ■ Vo)a;o, (6.9) 



dt 

we may find an equation for the time evolution of ^i. Using the fact the convective derivative 
is a derivation, we find 



vi X a;o + (V ■ xi)(vo x uq) + (V ■ cc;o)(xi x vq) + (V ■ vo)(u;o x xi). (6.10) 



-T^ = (Vi + (Xi ■ V)vo) XUJo + XiX {{ujo ■ V)vo - (V ■ Vo)t^o) 

dt 



In the second line the ordering of the symbols is meant only to indicate how the indices are 
wired up. The V must be understood to act to the right only on the velocity field vq. 
We now use the vector identity 

(x • a)(b X c) + (x ■ b)(c x a) + (x ■ c)(a x b) = x[a • (b x c)] (6-11) 

with X replaced by V (still acting only on vq) to find that 

^ = vi X u;o + V(xi ■ (vo X loq)) 
= Vi X cuo - V(vo ■ (Xi X Uo)) 
= Vi X cuo - V(vo ■ 6) 

= (VV^i) X cuo - (6 ■ V)vo (6.12) 



which is the same as (|5.15|) . (Again, in the first three hnes, V must be understood to act 
only on vq, even though it may be written to the left of other variables.) 

We can also check the consistency of the time evolution of the first order vorticity. From 
rSI) we find that 



so 



Ui = V X (Xi X Uq) 



V X I X CJo + V X Xi X 



dt \dt 7 V / 

It is not immediately obvious that (|6.14|) is compatible with the equation 

duji 



dt 



V X (Vo X CUi) + V X (Vi X UJo), 



which comes from applying d/dX to the vorticity evolution equation 

— = V X (v X cu). 
dt ^ ' 



The right hand sides of (|6.14|) and ( |6.15|) are equal only if 



9xi duji^ 

X CJo + Xi X — Vq X CJi — Vi X C(Jo, 



(6.13) 



(6.14) 



(6.15) 



(6.16) 



(6.17) 



dt " ^ dt 

is the gradient of something. Now by using ( |6.6| ), ( |6.13| ) and ( |6.15| ), we can write ( |6.17| ) as 

CUo X (V X (Xi X Vo)) + Xi X (V X (Vo X LJo)) + Vo X (V X (cJo X Xi)) 

-(t^o X Xi)(V ■ Vo) - (xi X Vo)(V ■ t^o) - (vo X a;o)(V • Xi). (6.18) 

Here we have added in a term (xi x vo)(V ■ tuo), which is of course identically zero, in order 
to preserve manifest cyclic symmetry of the terms. 

Now for any three vector fields a, b, c, we may verify that 

a X (V X (b X c)) + b X (V X (c X a)) + c X (V X (a X b)) 
-(a X b)(V ■ c) - (b X c)(V ■ a) - (c x a)(V ■ b) 

= V(a-(bxc)), (6.19) 



where V is acting on everything to its right. Applying this to ( |6.18D shows that it is a total 
derivative, and so the evolution equations are consistent. 
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VII. ILLUSTRATION 



As an illustration of the formalism consider waves propagating in the background flow 

f-y\ 



Vo 



X 



(7.1) 



V J 

This corresponds to the fluid rotating as a rigid body with angular frequency uJo/2. The 
perverse notation for the frequency arises because we have been using the symbol cu to denote 
vorticity, and V x vq = <^o z. 

To reduce notational clutter, in this section we will drop the suffix 1 from the fields 
and ^. It should still be borne in mind that they are first order quantities. 



A. Sound/Inertial- Wave Hybridization 



Our equations of motion are 



^ - (W X cuo) + {i ■ V)vo = 0. 



(7.2) 



They need to be supplemented with an initial condition that sets ^-uJo — 0. This orthogonal- 
ity is then preserved by the subsequent motion. Wc will ignore any effects due to gradients 
in po and (?. 

Take as an ansatz a plane wave in the frame rotating with the fluid. 



Here S^./ and are constant amplitudes. The primed unit vectors are 

ic' = xcos (^y^ i + ysin t, 
f = -xsin (^y) t + ycos (^y) ^' 



(7.3) 



(7.4) 
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and the primed coordinates 



X = X COS I — jt + ysm I — I t, 
y' = -X sin (^y ^ t + ycos (^y ) t, 



z = z. 



(7.5) 



The convective derivatives on tl) and on the components of ^ become 

= —iQi/j, 



dip 




H 


-m 






dt 


\ dt 



(7.6) 

For ^ itself we need to take note of the time dependence of the the unit vectors x', y', so we 
have 



dt ~ ' 



dt 



dt + y)^^' 



2^ . 

f ) ey)x' + {-^niy' + (^) ^.Oy' 



(7.7) 



Also we need 



K-V)v„ = -e,,(|l)x' + «„,(i|)y'. 



(7.8) 



The two off-diagonal uJo/2 terms add to get rid of the 1/2. The coupled equations 
therefore become 

/ —i^l —ujq —ikyiujQ \ I 'B.-j-i \ 
+ti;o —iVt +ikx'UJo Eyi = 0. (7.9) 
iky, [^-P)J \ ^ J 
For a solution to exist, the determinant of the matrix in ( [7.9|) must vanish. This gives 
the dispersion relation 



(^2_^2)/^_|^|2^ , ,.2,.2 , ,,2,.2 



ujQk^, + uj^ky, = 0, 
which for fixed k is a quadratic equation for Q^. 



(7.10) 
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Some insight into this dispersion relation can be obtained by letting (? 
limit the quadratic reduces to 



oo. In this 



(7.11) 



and so gives 



(7.12) 



This is the well-known dispersion relation for inertial waves in an incompressible fluid 20 



For these modes the restoring force comes entirely from angular momentum conservation. 
They are low frequency, ^ < Uq, oscillations and have a number of unusual features. In 
particular the frequency is independent of the magnitude of k, so the group velocity is 
perpendicular to the phase velocity — i.e. parallel to the wavecrests. At any particular 
frequency the disturbance spreads out from its source along a diabolic cone. 

The second root of the quadratic equation, ^ c^/c^ corresponds to conventional sound, 
and is lost to inflnity as becomes large. 

Now let us consider general values of c^. From the eigenmode equation we can solve for 
^ in terms of the amplitude of to get 

UJQ f -ky'Vt + ikx'UjQ 
- ^0 \ kx'Vl + iky^ujQ 



(7.13) 



This appears to be singular when VL^ approaches cUq, but, as we will see, this occurs only 
near k^' = ky/ = and the limit is smooth, the fluid rotating in circles in the x-y plane. 

From ^ we can flnd the velocity fleld, vi, and hence, by integration, the flrst order 
displacement fleld, xi, in the frame rotating with the background fluid^. We therefore flnd 
( ~kx'\ f —kx'UJo — iky/il\ 

—kytujQ + ikx'Vt 




kyl 



^0 
1]2 - ul 



\ 



J 



(7.14) 



^If xi = Xia;/x' + Xiy'y\ and vi = t^u-'x' + viy'y\ then ( [6.6D reduces to Vix'^y' = f—^^) 
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It is now straightforward to verify that we recover ^ from ^ = (xi x uo). We also verify 
that the correction to potential flow is 0{uJo/Q) when fl ^ uq. 



B. Poincare Waves 

If we restrict ourselves waves with kz = 0, then setting the determinant to zero gives 

fi2 



C2 



(u'^ + c^k'^ -n'^) =0. (7.15) 



We therefore have two classes of modes: those with zero frequency, and those with a gapped 
dispersion relation 

n' = u;l + c\kl, + kl). (7.16) 



The former are z independent (Taylor- column pT| , p2| ) geostrophic flows where pressure 
gradients are balanced against a Coriolis force. The gapped modes are the Poincare waves 



We can obtain the Poincare modes by considering the motion directly in the x', y' frame. 
The effect of the frame rotation produces a Coriolis force, and so the equation of motion is 

dvix' 

—g^ = i^OVly' - C Vx'Pl, 

^'"^y' 2V7 /^7 1 ^7^ 

= -ujqVix' - c Vy'Pi. (7.17) 
To solve we need to combine this with continuity equation 

^ + PO {Vx'Vlx' + Vy'Vly') = 0. (7.18) 

For waves travelling in the x' direction we flnd 



vix' = Acos{kx' — Qt) 

UJc 

oo 



viyi = A (^"^^ sm{kx' — Qt), 
p^= A(^]cos{kx' -nt), (7.19) 
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together with the dispersion relation ( [7.16| ). 

We also find the displacements of the particles to be 



xix' = —A ( — ) sm{kx' — Qt), 



x^y, = A (^^j cos(A;x' - nt). (7.20) 

When k = 0, we have Q = uq, and the particles move in circles in the x, y plane. This 
limiting motion is solenoidal and coincides with the kx' = ky/ = limit of the incompressible 
fiuid inertial waves. 

We may now make contact with our ip, ^ formalism by writing 

\ / x',y' 

where 

fc^k\ 

^ = Al — \sm{kx' -Qt). (7.22) 

The relation ,^ = xi x wq, which was derived in the x, y inertial frame, continues to hold 
in the rotating frame without modification. So we have 

^ = xi X tuo = uJoz{xiy>±' - xix'-y'). (7.23) 
We can combine our expression for ip with this to get 

^ + ^ I cos(A;a;' - nt) 
v^y, = +^y, = A (^^^ sin(A;x - nt). (7.24) 

Since the factor in parenthesis in the first line is seen to be unity by use of the dispersion 
relation, we recover the earlier expression for vi, and confirm that the gauge invariant 
decomposition works as advertised. Again we see that the velocity field ^, which arises from 
angular momentum conservation, is smaller than the pressure induced fiow, Vip, by a factor 

of UJq/Q. 
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VIII. DISCUSSION 



The central idea in this paper is the decomposition Vi = V^' + C of the velocity pertur- 
bation into a potential flow and a correction required by angular momentum conservation. 
This decomposition is motivated by the Clebsch formalism, but does not depend on it. From 
the decomposition we see that corrections to the acoustic metric equation depend only on 
the ratio of the frequency of the sound wave to the frequency, cuq/S, of the background fluid 
rotation. This rotation frequency is determined by the antisymmetric part, {diVj — djVi)/2, 
of the velocity inhomogeneity. The symmetric part, the rate of strain — {diVj + djVi)/2, 
can be large and the correction remain small. This is not unreasonable because the acoustic 
metric equation is exact for any potential background flow — no matter how inhomogeneous. 

At low frequencies the correction ^ = Xi x a;o ceases to be negligible. In this regime the 
sound waves hybridize with whichever of the many other modes available to a fluid with 
vorticity happen to have comparable frequency. The hybridization may lead to a spectral 
gap, as with the Poincare waves, to birefringence, and to other phenomena which show that 
the acoustic metric is no longer all that is needed to describe sound propagation. 
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